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The Global Positioning System enables positioning in 3 dimensions about our planet. The 
signals sent by twenty-four satellites travel through our atmosphere down to the small receivers 
used by the civilian community and the military. These signals are electromagnetic waves; they 
do not travel in a straight line because of varying meteorological conditions. Instead, they bend 
towards the surface and travel slower than in a vacuum. 

A small GPS receiver was launched in the satellite GPS/MET in 1995. It became possible 
to perform radio occultalions around the Earth: the source - one of the 24 GPS satellites - is 
seen by the receiver as il rises or sets around the other side of the Earth. W hen the source 
disappears, the receiver progressively loses the signals. The delay induced by the atmosphere 
can be infered by measuring measuring accurately the time delay between the emission and 
the reception of the sign? 1. 

In a previous study we had shown one method to combine GPS/MET data with numerical 
weather forecasts in a statiscally optimal way to yield improved profiles of temperature. 

We use these observations and assimilate them in the Data Assimilation Office Finite 
Volume Data Assimilation System (FVDAS). The information brought by the temperature 
profiles is used by the FVDAS to perform global three dimensional weather analyses of the 
atmosphere. We evaluate the impact of the GPS/MET data on the global analyses and the 
weather forecasts. 
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1 Introduction 

The GPS/MET 1195 experiment demonstrated the feasability of remote 
soundings of the Earth’s atmosphere using the GPS radio occultation tech- 
nique [5]. Subsequent studies assessed the information content of such mea- 
surements (e.g. [3]) Several groups have attempted to assimilate the GPS ra- 
dio occultation measurements into a global Data Assimilation System (DAS). 
These attempts followed in general the methodology proposed by Eyre [2], 
and range from direct four-dimensional variational assimilation of bending 
angles to assimilation of inverted profiles of temperature and/or humidity. 
We present here an hybrid approach. We attempted to take advantage of the 
accurracy of curren . w r eather models to constrain one-dimensional variational 
analyses using GPS refractivity. Retrieved profiles of atmospheric parameters 
were then assimilated like other types of observations in a global DAS for is- 
suing the next forecast period. 


2 One-dimensional Variational (1DVAR) Analyses with 
Refractivity 

The application of one-dimensional variational analysis (ID VAR) to the re- 
trieval of geophysh al parameters from GPS radio occultation has already 
been discussed by various authors (e.g. [3], [9]). We briefly describe the main 
ideas. Technical details can be found in [10]. 

2.1 1DVAR Analysis Theory 

Air refractivity N for GPS signals frequencies (L-band) can be written (ne- 
glecting scattering and ionospheric contribution) 

N = b,^ + b 2 ^ 


( 1 ) 
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[11], where P is the air pressure in hPa, T the temperature in K, P w 
the water vapor partial pressure in hPa, b\ = 77.6 N-unit-K-hPa -1 , and 
b 2 = 3.73xl0 5 N-unit K 2 hPa -1 . It is obvious from equation (1) that the 
direct inversion of refractivity measurement into temperature, humidity and 
pressure is not possible without assumptions. One way to address the tem- 
perature/humidity ambiguity is to use physical constraints via a 1DVAR 
approach. A priori information (or background) is combined with GPS re- 
fractivity observations in a statistically optimal way to yield analyzed profiles 
and temperature and humidity. An important hypothesis is that all errors are 
supposed to be non-biased and normally distributed. 

2.2 GPS/MET Observations 

We used 797 occupations collected during ‘Prime Time 2’ of the GPS/MET 
1995 experiment, when the GPS Anti-Spoofing encryption had been turned 
off from June 21st until July 4th, 1995. This, along with changes in the 
software of the receiver enabled tracking of the GPS occupations deep in 
the troposphere. In order to use these obervations in the 1DVAR scheme 
we assumed the following standard deviation error in refractivity: 1% below 
5 km aPitude, 0.2% between 5 km and 30 km altitude. We discarded the 
refractivity observations above 30 km aPitude. 

2.3 Background Information 

The Data Assimilation Office (DAO) Finite Volume Data Assimilation Sys- 
tem (FVDAS) [1] provided the 6-hour forecasts used as a background in 
the ID VAR implementation. The forecasts were bilinearly interpolated from 
2x2.5 degrees grid to the location of the GPS occupation. No time interpo- 
lation was performed. The forecast fields (from the surface to the 0.01 hPa 
level) were interpolated to 46 levels for the ID VAR analysis. 

2.4 GPS 1DVAR Results 

Due to a positive bias in the background minus observation refractivity the 
analyzed 1DVAR humidity profiles presented systematic departure from the 
original 6-hour forecast used to initiate the retrieval. These departures were 
inconsistent with humidity measurements reported by nearby radiosondes 
[10]. We hence decided to investigate this subject further before assimilating 
these analyzed humidity profiles into the FVDAS. 

Figure 1 shows temperature comparisons with nearby RS (less than 
280 km, +/- 3 hours). The 1DVAR analysis has smaller bias than the back- 
ground above the 200 hPa level. The standard deviation of the background is 
larger than the ID VAR analysis between 850 hPa and 20 hPa. This suggests 
that the 1DVAR analyzed temperature profiles may provide useful informa- 
tion to a global Data Assimilation System (DAS). 
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Fig. 1 . Temperatun comparisons of 6-hour forecasts and GPS 1DVAR with nearby 
radiosondes; numbej of comparisons on the right-hand side 


3 Assimilation of GPS 1DVAR Profiles in the Finite 
Volume Data Assimilation System (FVDAS) 

3.1 Description of the FVDAS and Link with GPS 1DVAR 
Analysis 

The FVDAS relies on a NASA/NCAR (National Aeronautics and Space Ad- 
ministration/National Center for Atmospheric Research) Finite Volume Gen- 
eral Circulation Model (GCM). The GCM dynamical core uses a Lagrangian 
vertical coordinate system [7]. The assimilation part of the FVDAS is a 
Physical-space Statistical Analysis System (PSAS) [1]. The FVDAS assimi- 
lates conventional observations such as radiosondes (RS) as well as 1DVAR 
cloud-cleared retrieved profiles of temperature and humidity derived from 
TIROS Operational Vertical Sounder (TOVS) level-lb brightness tempera- 
tures [4]. 

When running he FVDAS from June 21st, 1995 to July 4th, 1995, the 
latest available forecast was used to perform the GPS 1DVAR analysis. Then, 
the GPS ID VAR analyzed temperature profiles were assimilated in the FV- 
DAS using PSAS. They were assigned the same error characteristics as ra- 
diosondes. Due to eignificant refractivity differences between 6-hour forecast 
background and observation [10], we chose to restrict the assimilation of 
1DVAR retrievals to the region 500 hPa-30 hPa in the [30°S,30°N] latitude 
band. At other latitudes, 1DVAR retrievals were assimilated between the 
surface and the 30 liPa level. We refer to this experiment as 'GPS 5 . 

In order to evaluate the impact of the GPS 1DVAR analyses on the FV- 
DAS we also performed a ‘CONTROL’ experiment with no GPS data. 
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3.2 Impact on Analyses 

Without any other observation other than those assimilated in the two ex- 
periments, it is difficult to determine which analysis is the closest to the 
atmospheric “truth”. Figure 2 shows the difference between an analysis with 
GPS data minus the corresponding ‘CONTROL 5 analysis for 1995-06-29 
at 00GMT. In this example we used the same first guess for perform both 
analyses. Consequently, the differences we see are solely due to the GPS tem- 
perature observations assimilated at 00GMT on 1995-06-29. The differences 
in the analyzed winds enables to assess how the GPS temperature observa- 
tions changed the atmospheric circulation through the mass-wind balance. 
For example, negative temperature increments corresponding to negative 
geopotential heights increments induced cyclonic- type increments in winds, 
i.e. horizontal circulation differences oriented clockwise in the Southern hemi- 
sphere. 



Fig. 2. Difference ‘GPS’ analysis minus ‘CONTROL’ at 250 hPa. Contours desig- 
nate geopotentiai heights, arrows represent wind. 


3.3 Impact on 6-hour Forecasts 

We compared the 6-hour forecasts from the two experiments with radiosonde 
observations. These observations may be considered as an independent vali- 
dation. Figure 3 shows zonal averages of the standard deviation of the obser- 
vation minus forecast for the meridional wind. The ‘CONTROL’ has larger 
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errors in the Southern hemisphere than in the Northern hemisphere where it 
is more constrained by in situ observations. The shaded areas in the ‘GPS’ 
minus ‘CONTROL’ indicate reduced errors in the ‘GPS’ 6-hour forecast as 
compared to the ‘C ONTROL’. Differences are more apparent in the Southern 
hemisphere. 


CONTROL 



CONTROL-GPS 



Fig. 3. Radiosonde observations minus 6-hour forecast standard deviation error for 
meridional wind. Negative coutours (dotted lines) in the ‘CONTROL’ minus ‘GPS’ 
suggest a negative impact of the GPS data; Shaded areas indicate a slight positive 
impact. 


3.4 Impact on 5 day Forecasts 

We ran a series of l -day forecasts from daily 00 GMT analyses from June 21st 
until July 4th, 1995. Forecast scores were computed in the form of Root Mean 
Square (RMS) and anomaly correlation for sea level pressure and SOOhPa 
geopotential height . We used three different verifying analyses: (a) own anal- 
ysis for each experiment, (b) ‘CONTROL’ analyses, and (c) European Centre 
for Medium-range Weather Forecasts (ECMWF) operational 1995 analyses. 
The GPS data induced slight differences on forecast scores of individual fore- 
casts. The average of the 14 forecasts showed very little difference between 
‘CONTROL’ and ‘GPS’. The small number of forecasts issued and the lim- 
ited amount of GFS data assimilated does not enable us to fully assess the 
impact of GPS refi activity data in a global data assimilation system. 


4 Conclusions and Future Work 

We have implemenl ed a ID VAR scheme to analyze temperature and humidity 
profiles using GPS radio occultation refractivity. This scheme was incorpo- 
rated in the assimi ation cycle of the FVDAS and we performed assimilation 
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experiments using GPS/MET refractivity. We found mixed results in terms 
of impact of the GPS data on forecasts. Clearly, the limited amount of GPS 
observations l im its the significance of the present work and more GPS data 
are needed to assess their impact on a global data assimilation system. 

Even in a 1DVAR framework with limiting hypotheses such as no account 
of the atmospheric horizontal structure, we found a positive impact of GPS 
refractivity on the analysis of temperature. However, the humidity retrievals 
were not satisfactory partly because of a “refractivity bias”. The source of 
this bias observed also in the Champ radio occultation measurements [8] 
has to be understood, because it plays even a more critical role in a direct 
assimilation where analysis increments cannot be separated. Similar assimila- 
tion experiments will be performed using the Champ and SAC-C GPS radio 
occultation observations, and the refractivity bias will be investigated. 
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